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Abstract 

The  interaction  between  a  shock  wave  (attached  to  a  wedge)  and  small  amplitude,  three- 
dimensional  disturbances  of  a  uniform,  supersonic,  freestream  flow  are  investigated.  The 
paper  extends  the  two-dimensional  study  of  Duck  et  aP,  through  the  use  of  vector  potentials, 
which  render  the  problem  tractable  by  the  same  techniques  as  in  the  two-dimensional  case,  in 
particular  by  expansion  of  the  solution  by  means  of  a  Fourier-Bessel  series,  in  appropriately 
chosen  coordinates. 

Results  are  presented  for  specific  classes  of  freestream  disturbances,  and  the  study  shows 
conclusively  that  the  shock  is  stable  to  all  classes  of  disturbances  (i.e.  time  periodic  per¬ 
turbations  to  the  shock  do  not  grow  downstream),  provided  the  flow  downstream  of  the 
shock  is  supersonic  (loosely  corresponding  to  the  weak  shock  solution).  This  is  shown  from 
our  numerical  results  and  also  by  asymptotic  analysis  of  the  Fourier-Bessel  series,  valid  far 
downstream  of  the  shock. 
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1  Introduction 

The  interaction  of  freestream  disturbances  with  shock  waves  is  an  important  practical  prob¬ 
lem,  with  particular  regard  to  the  operation  of  high-speed  flight  vehicles.  Of  particular 
interest  are  the  consequences  of  these  disturbances  on  boundary-layer  receptivity,  and  this 
paper  may  be  regarded  as  a  step  in  this  process. 

In  the  case  of  an  isolated  shock,  Moore  Ribner  ^  and  McKenzie  &:  Westphal  ^  showed 
that  a  single  wavelength  entropy,  acoustic  or  vorticity  wave,  upstream  of  the  shock,  produces 
a  disturbance  comprising  a  mixture  of  all  three  classes  downstream  of  the  shock.  When  the 
shock  is  associated  with  the  flow  past  a  rigid  body,  the  situation  is  much  more  complicated, 
due  to  the  reflection  and  refraction  of  disturbances  between  the  body  surface  and  the  shock 
itself,  the  latter  distorting  as  a  result,  a  process  which  further  complicates  the  flow  structure. 
This  process  was  considered  by  Carrier  ®  and  Van  Dyke  ®  with  particular  regard  to  the 
problem  of  supersonic  flow  past  a  wedge  performing  small  amplitude  oscillations,  the  shock 
remaining  attached  to  the  wedge  tip.  These  problems  raise  questions  regarding  the  stability 
of  the  shock,  and  this  aspect  has  been  considered  in  the  two-dimensional  context  (associated 
with  wedge  flows,  the  shock  remaining  attached  to  the  wedge  tip  at  all  times)  by  Levinson’^, 
Carrier  Henderson  &  Atkinson  Rusanov  &  Sharakshannae  and  Salas  &;  Morgan^^  The 
overall  conclusion  is  that  if  the  flow  downstream  of  the  shock  is  subsonic  (loosely  classified 
as  the  strong  shock  solution),  then  the  shock  is  unstable,  in  so  far  as  disturbances  grow 
downstream.  If,  on  the  other  hand,  the  flow  behind  the  shock  is  supersonic  (loosely  classified 
as  the  weak  shock  solution),  then  the  shock  is  stable,  with  disturbances  not  growing  in 
amplitude  downstream. 

This  aspect  was  considered  in  some  detail  by  Duck  et  al.  again  for  the  case  of  super¬ 
sonic  flow  past  wedges,  subject  to  two-dimensional  disturbances.  This  latter  paper  (where 
there  is  also  a  more  detailed  summary  of  previous  work  in  this  area)  showed  that  two- 
dimensional  disturbances,  impinging  upon  the  weak  shock  solution,  do  not  grow  downstream, 
but  rather  produce  constant  amplitude  or  weakly  decaying  oscillatory  waves  downstream. 
The  conditions  which  distinguish  between  the  two  types  of  behaviours  were  determined  and 
a  decomposition  of  the  far-field,  suggested  by  the  analytic  solution,  was  made.  The  result 
was  that  there  is  a  shadow  region  dependent  on  the  incoming  disturbance.  If  the  direc¬ 
tion  of  propagation  of  the  plane-wave  acoustic  disturbance  generated  at  the  shock  by  the 
plane-wave  incoming  disturbance  intersects  the  wedge  surface,  then  there  is  a  non-decaying 
oscillatory  pressure  disturbance  on  the  wedge  surface.  In  any  case,  there  are  two  weakly 
decaying  single-wavelength  oscillatory  disturbances  on  the  wedge  surface  that  are  due  to 
the  requirement  that  the  shock  remain  attached.  One  is  traveling  at  the  mean-flow  speed 
plus  the  speed  of  sound,  and  the  other  travels  at  the  mean-flow  speed  minus  the  speed  of 
sound.  Thus,  these  are  clearly  regular  acoustic  waves  which  emanate  from  a  point  source 
(the  wedge  apex)  of  fixed  temporal  frequency  as  prescribed  by  the  incoming  disturbance.  In 
addition  to  these  more  obvious  effects,  there  is  also  a  local  field  effect  which  is  broad  based 
in  wavelength.  While  the  first  three  modes  of  pressure  disturbances  at  the  wedge  surface 
can  indeed  interact  with  the  boundary  layer,  it  is  the  last  disturbance  with  the  inclusion  of 
short  wavelength  disturbances  that  should  have  the  most  effect  on  the  growing  boundary 
layer.  In  this  latter  paper,  the  side  issue  of  shock  stability  to  two-dimensional  disturbances 
was  also  addressed.  The  aim  of  the  present  paper  is  to  extend  the  ideas  of  Duck  et  al.  ^  to 
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three-dimensional  disturbances;  indeed,  little  three-dimensional  work  of  this  type  appears  to 
have  been  undertaken  in  the  past. 


2  Formulation 

Throughout  this  paper,  we  use  subscript  1  to  denote  quantities  upstream  of  the  shock,  and 
subscript  2  to  denote  downstream  quantities.  The  wedge  is  taken  to  make  an  angle  9  with 
respect  to  the  oncoming  flow,  with  {x*,y*,z*)  coordinates  parallel  and  perpendicular  to  the 
upstream  flow,  with  z*  being  the  “crossflow”  direction.  The  upstream  basic  flow  velocity 
has  magnitude  f/*,  with  Mach  number  Mi  and  density  and  temperature  Tj*.  The  ratio 
of  specific  heats  7  is  assumed  to  be  constant,  and  throughout  this  paper  all  our  numerical 
results  take  7  =  1.4.  The  velocity  vector  is  written  as  t/i*(u,  u,  re),  with  respect  to  (a;*,  y*,  z*) 
coordinates,  and  the  density  is  then  written  p*p,  pressure  as  p\R*TiP  (where  R*  denotes  the 
gas  constant)  and  the  temperature  is  T^T. 

The  Rankine-Hugoniot  relations  provide  a  link  between  conditions  upstream  and  down¬ 
stream  of  the  shock.  For  the  basic  (steady)  flow,  when  the  downstream  flow  is  uniform  and 
parallel  to  the  wedge  surface,  the  following  classical  result  is  obtained  (see,  for  example, 
Liepmann  and  Roshko  ^^): 

tan(/?  -  ^)  _  il2  _  pi  _  (7  -  l)Ml  sin^  /?  -f  2 

tany0  ui~  P2~  (7-f  l)M2sin2/?  ’  ^ 

where  (3  is  the  angle  between  the  shock  and  the  wedge  centerline  (i.e.,  y*  —  0)  and  ui  and 
U2  are  the  non-dimensional  velocity  components  perpendicular  to  the  shock. 

The  relationship  (2.1)  yields  two  possible  values  for  /?,  for  a  given  value  of  6  (for  9  < 
^max(Mi))  with  the  so-called  “weak  shock”  solution  generally  being  characterized  by  super¬ 
sonic  flow  downstream  of  the  shock,  whilst  the  other  “strong  shock”  solution  is  generally 
characterized  by  downstream  subsonic  flow.  Note,  however  that  there  does  exist  a  small, 
weak  shock  regime,  close  to  ^  =  9max  where  the  downstream  flow  is  subsonic,  i.e.  the 
downstream  sonic  line  does  not  quite  coincide  with  9max- 

We  shall  be  concerned  with  the  effect  of  small  amplitude  disturbances,  which  are  intro¬ 
duced  into  the  flow  ahead  of  the  shock.  The  corresponding  two-dimensional  results,  as  stud¬ 
ied  by  McKenzie  and  Westphal  for  example,  may  be  extended  into  the  three-dimensional 
contact,  in  quite  a  straightforward  manner,  when  the  z*  variation  is  taken  to  be  periodic  (al¬ 
though  other,  more  general  z*  variations  could  be  accomplished  using  Fourier  Transforms). 
We  take  e  (<C  1)  to  be  a  measure  of  the  amplitude  of  the  freestream  disturbance,  and  it 
then  turns  out  that  the  upstream  disturbance  may  be  classified  into  three  distinct  classes, 
just  as  in  the  two-dimensional  case. 

(i)  Acoustic  waves:  these  are  characterized  by  having  a  pressure  perturbation,  with 
corresponding  perturbations  in  velocity,  density  and  temperature,  but  no  change  in  entropy 
or  vorticity,  to  0(e).  Disturbances  of  this  class  take  the  form 

p  =  1 -f- ^ -f- 0(6^),  (2.2) 
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u  =  1  — 


aisE 


V  = 


w  — 


'yMf{ai  +  u) 
a2eE 


+  0(6% 


'yM^{ai  +  w) 

A 

a^eE 

+  w) 


+  0{e% 
+  0{e% 


and 


p  =  1  — -  +  O(e^), 
7 


T  =  1  -  iil— ^  +  0(£2), 
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where  E  is  the  normal  mode  exponential,  i.e. 

E  =  exp{miXi  +  ^0^22/1  +  io^sz  +  icot}y 

with 

w  =  -a,  ±  -^[aj  +  o|  +  41'/* 


(2.3) 

(2.4) 

(2.5) 

(2.6) 

(2.7) 

(2.8) 

(2.9) 


being  the  frequency  for  the  given  wavenumbers.  In  the  above,  x\  and  j/i  are  parallel 
and  perpendicular  to  the  upstream  flow  respectively  (both  perpendicular  to  2);  the  non- 
dimensionalisation  is  carried  out  using  one  of  the  wavelengths  of  the  disturbances  as  the 
typical  lengthscale,  e.g.  by  setting  ai  to  unity.  The  modes  with  the  positive  sign  above 
in  (2.9)  are  usually  referred  to  as  the  slow  modes,  whilst  those  with  the  negative  sign  are 
usually  referred  to  as  the  fast  modes. 

(ii)  Vorti city  waves:  these  are  characterized  by  having  no  density,  temperature  or  pressure 
disturbances  to  0{e),  and  thus  the  upstream  flow  takes  the  form 


u  =  1  +  eO. 


V  =  e 


-n 


ai 


al  +  ol\ 
w  =  eilx- 


+  0(s% 

(2.10) 

aj  -t 

rt  ""  1 

E  +  0{e^)^ 

(2.11) 

^E  +  0(e^), 

2  ^  ^3 

(2.12) 

—  1  +  0{e^), 

(2.13) 

}  =  —ai. 

(2.14) 

and 

with 


The  exponential  term  E  is  unaltered  from  (2.8),  and  the  parameters  and  Clx  must  be 
specified.  In  comparison  with  the  two-dimensional  case,  the  three-dimensional  case  allows 
for  an  extra  vorticity  mode. 
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(iii)  Entropy  waves:  these  have  no  disturbance  pressure  or  velocity  components  to  0(e), 
and  may  therefore  be  written  in  the  form 

p  —  \  eE  O(e^),  (2.15) 

T  =  l-  eE  +  0{e^),  (2.16) 

u,p  =  l  +  0(e'‘),  (2.17) 

and 

v,w  =  0{£^),  (2.18) 

with  uj  defined  by  (2.14)  and  E  by  (2.8),  again. 

In  the  following  section  of  the  paper,  we  consider  the  effects  that  these  waves  have  on  a 
shock  wave  attached  to  a  wedge.  The  complication  introduced  by  the  presence  of  the  wedge 
is  that  a  single  mode  of  the  above  type,  upstream  of  the  shock,  will  trigger  modes  of  all  three 
types,  of  all  wavenumbers  behind  the  shock;  the  only  quantities  preserved  across  the  shock 
(provided  e  ^  1)  are  the  frequency  parameter  (u  and  the  spanwise  wavenumber  0:3. 


3  The  Downstream  Solution 

The  flow  downstream  of  the  shock  comprises  the  uniform  flow  solution  plus  a  small  ampli¬ 
tude  (0(e))  perturbation,  triggered  by  the  freestream  disturbances  described  in  the  previous 
section.  We  take  non-dimensional  coordinates  parallel  and  perpendicular  to  the  wedge  (per¬ 
pendicular  to  the  direction  in  both  cases)  as  X2  and  y2i  respectively,  with  corresponding 
velocity  components  (u2,V2). 

We  now  write  the  solution  downstream  of  the  shock  in  the  form 


U2  —  U2  -|-  EU  0(e^), 

(3.1) 

V2  =  EV  +  0(e^), 

(3.2) 

W  =  EW  +  0(e^), 

(3.3) 

P  —  P2  £p  0(e'^), 

(3.4) 

p-p2+Ep  +  0(e^), 

(3.5) 

T  =  T2  +  Ef  +  0(E^). 

(3.6) 

Substitution  of  these  expansions  into  the  governing  (inviscid)  equations  of  motion,  continuity 
and  energy  equations  and  equations  of  state  then  leads  to 


Pt  +  U2PX2  +  P2UX2  +  P2^y2  H"  P2Wz  —  0, 


(3.7) 


P2{ut  +  CiUij}  +  —  0, 


(3.8) 
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(3.9) 


and 


P2{vt  +  U2VX2}  +  -  Oj 

P2{wt  +  U2W^J  +  -4t2Pz  =  0, 

im 


P2{Tt  +  U2TX2} 
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{Pt  +  U2PX2}  =  0, 


(3.10) 

(3.11) 


p  =  p2f  +  T2P. 


(3.12) 


Much  of  the  success  of  the  approach  used  by  Carrier  Van  Dyke  ®  and  indeed  by  Duck 
et  al.  ^  was  due  to  the  ability  of  being  able  to  split  the  solution  into  two  components,  the  one 
including  a  velocity  potential  (which  represented  the  acoustic  wave  component  of  the  flow), 
and  the  other  a  stream  function  (which  represented  the  vorticity  wave  component  of  the 
flow).  It  is  clearly  not  possible  to  use  the  latter  in  the  context  of  three-dimensional  flows, 
however  an  alternative  is  the  introduction  of  a  vector  potential.  This  concept,  originally 
due  to  Poincare  has  been  used  in  a  number  of  fluid  mechanics  investigations  over  the 
years,  although  the  appropriate  boundary  conditions  have  been  the  subject  of  discussion,  as 
detailed  by  Aziz  and  Heliums 

Specifically  we  write  the  perturbation  velocity  vector  u  =  (u,v,w)  in  the  form 


u  =  Vp  -t-  V  A  E, 


(3.13) 


where  4>  is  the  velocity  potential,  and  E  the  vector  potential.  The  above  represents  a  non¬ 
unique  representation  for  u,  and  as  such  an  additional  relationship  may  be  (arbitrarily) 
specified.  The  most  popular  choice,  and  the  one  that  we  pursue  here,  is  that  the  vector 
potential  can  be  required  to  be  solenoidal,  i.e. 


<1 

II 

0 

(3.14) 

This  leads  to  the  vorticity  vector  then  being  merely  the  quantity  — V^E. 
(3.13)  into  (3.8)  -  (3.10)  leads  to  the  equations 

Substitution  of 

(3.15) 

and 

W^Et  +  U2E,,}  =  0. 

(3.16) 

It  then  follows  that  p  satisfies 

(3.17) 

and  elimination  of  the  dependent  variables  p  and  T  gives  the  equation  for  </> 

=  —{(ptt  +  2C/20a;t  +  U24>X2X2\i 

CZ-- 

(3.18) 
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where 


(3.19) 


If  we  write  E  =  then  on  account  of  our  comments  above,  we  impose  the 

condition 

Ei^J  +  Ell'>  +  Ei^^  =  0.  (3.20) 

A  further  quantity  must  also  be  introduced,  namely  'tl)'{y2,z,t),  which  describes  the  displace¬ 
ment  of  the  shock  wave  from  its  undisturbed  state. 

This  type  of  formulation  then  allows  us  to  follow  closely  the  two-dimensional  approach 
adopted  by  Duck  et  al.  h  It  is  now  possible  to  write  the  general  solution  for  <^,  and  the 
ew  in  the  following  form,  assuming  boundedness  at  the  apex  of  the  wedge: 


.  .  iuUoxo  ,  .  OO 

twt-  >  ^  I  +tasz  _ 


^  cosh(r/02)  +  K  smh{v92)}J,y{ker), 


(3.21) 


.  ,  tu/xo  tojAyoa?  ,  .  CX) 


(3.22) 


f  =  e 


y  ^  diiJv{ke^y2^  1 


(3.23) 


where  we  have  written 


(3.24) 


A  (jJ 

Ul  -  al 


(3.25) 


(3.26) 


kl  =k^  +  ^, 
^2’ 

A  =  cot(/3  —  6), 

A  j  A 

e  = 

r'^  =  x\-  PY2, 


(3.27) 

(3.28) 

(3.29) 

(3.30) 


tanh  $2  —  i^y2lx2- 

If  we  impose  impermeability  on  the  wedge  surface,  we  immediately  require 


(3.31) 


bu  =  0  Vi/, 


(3.32) 
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whilst  if  the  shock  remains  attached  at  the  tip,  then  do  =  0.  As  noted  by  Duck  et  al.  '*  (for 
example),  it  is  also  reasonable  to  set  the  z/  =  0  coefficients  of  the  other  terms  in  the  series 
to  zero,  which  correspond  to  transient-type  modes,  and  so  we  write 


ao  =  c^'>  =  =  cP  =  0. 

The  key  results  for  the  velocity  components  and  pressure  are  then 


(3.33) 


“  L  r  2im 


«  =  -e 


u=0  2  [keP'^di 


%i,  +  ay_i(l  -f  5;/-i,o)  —  «£/+!  cosh  v'02Ju{ker) 


rL  f  2fAd;a 


^oL2  V  kM2f3^ 


+  e  (eg,  -  c2i(l  -b  <J,_i,o))  -  f/342) 


1  (a^_i -t- a,,+i)sinh(i/02)d'./(fcer) 

z/=0  ^ 


V  =  e  “5 


.  .  iwx2  twAya,  ,  .  OO 


E  +  Mkeiy2), 

2=n  L  U2  J 


.  iUoa/zEo  ,  •  OO 

+»‘^3Z 

w  —  e  2  > 


y~)  i^ai,J^{kf.f’)  cosh 


(3.34) 


(3.35) 


•  .  ,  •  tazreo  twAt/oa?  OO 

«u;*+ta3j5 - - A'  ■■ 


'■  S +  T  +  <5.-..o)))  I  UkM, 


S  C/2 


(3.36) 


~  '1P2ke  iu>t-'-§¥^+ia:,z  \  2iCb  ^  ^  Uo 

P  ~  ~y, - ®  ^  *  X)  - “i4-l(l  +  ^,/-l,o) - fliz+l  COsh(l/02)«/./(^e>')) 

^,-0  [fceP  J 

(3.37) 

where  c_i  =  cij  =  d_i  =  0,  and  where  Sn,m  is  the  Kronecker  delta. 

Although  the  solutions  above  satisfy  the  impermeability  and  apex  conditions,  as  yet 
these  solutions  do  not  satisfy  the  unsteady  Rankine-Hugoniot  conditions  (McKenzie  and 
Westphal'*)  on  the  shock;  following  (and  extending)  the  results  of  Duck  et  al.  \  these  take 
the  form 


p  zr - e  ^2  ""s 

2a. 


u'  -b  U2—  =  ( 1 


M  nz  ,  1 


R'  -I — {u,pi  -b  «i), 

P‘2 


/  1_  P2  .  a*  P2  ,  1 _  Pi 

^2  ”b  77^2 - 1”  - - —  U.  -| - UiPi  H - o~~  1 

2  p2  27U2P2  2  2'yMiUi 


(3.38) 


(3.39) 


<  P2 


- ——  _| - - 

(7-l)u2p2  (7-1)^^2P2 


a;  P2 


(3.40) 


■(7-l)M2u2’*’(7-l)Mfu2’ 
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and 


(3.41) 


v'^  =  v[-{-  (ui  -  'U2)V’^, 

w'2  =  w[  +  (ui  -  ■u2)V’z,  (3.42) 

with 

i?' =  (3.43) 

In  equations  (3.38)-(3.43),  a  denotes  the  unperturbed  distance  along  the  shock,  {u',v',w') 
denote  the  velocity  perturbations  perpendicular  to  the  undisturbed  shock,  parallel  to  the 
shock  and  perpendicular  to  the  2;  direction,  and  parallel  to  the  z  direction,  respectively,  and 
(u,u,0)  denote  the  corresponding  base  flow  velocity  components.  Eliminating  p2  from  the 
above  equations,  and  then  utilizing  (3.34)-(3.37),  yields  the  following  four  equations  for  each 
of  the  u  >0: 


cos{f3 


2iu)U2 


du  +  <*i/-l(l  +  ^i/-l,o)  ~ 


cosh  pOq 


+ 


+ 1  (cS.  -  41(1  + 

.  A  A 


+  sm{0  -  0)  |-|-(a^_i  +  a,.+i)  sinh(z/0o)  + 


-(«!  -  U2)  sm{0  -  0)j 


dp  i  (^i/+i  ““  (3.44) 


4li 


+ 


j^sin(/?  -  9) 

A  . 

ke  I  2iXcba 


k^ 

T 


2iuU2 


dp  +  dv-i{l  +  Si,  -1,0)  —  Olz+l 


2  V  hU20^ 
h0 


ke0'^d 

4"^  +  i  (4?i  -  cSi(i  +  <^.-1,0)))  -  i0c?^ 


cosh  pOq 


—  cos(0  —  0) 


{a„^i  +  a„+i)  sinh(i/0o)  +  +  77-4^^ 

Z  U2 


kp 


2iu}U2X 


ke/S'^ds 


dnd  +  i  (d^+1  ~  d^_i(l  +  (J;,_i^o)) 


+  iujdv 


+ 


{^111 


'2iu  U2  .  U2  \ 

-j—^d„  H - +  (),._i,o) - a^+i  >  cosh(i/6'o) 

kep^  <^s  Cls  )  ^ 


4^) 


(3.45) 


and 


ias  cosh(?/^o)av— ^4^^ f  4^^+|  (cli  -  cj,l\)(l  +  4-i,o))-ia3(«i-U2)d,,  = 


U2  2  ^  fceC/2/3' 


(3.46) 
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These  equations  must  also  be  augmented  with  the  condition  arising  from  the  solenoidal 
condition  on  the  vector  potential  namely 


iu)  ke ,  2iXi!bas 
+  ir(' 


U2'''  '2'  k,U2^^ 

Terms  used  in  the  above  are 


4^^  +  i  (4+1  -  4-1  )(i  +  <^i^-i,o))  +  *"34^^  =  0- 


tanh  Oo  =  P  tan(/3  —  9) , 

'^■4' 

'4^=4i-4. 


Bi  =  l  + 


2ui' 


Bo  = 


U\ 


1 


L(7-1)“2  “2J 

■yds 


1--  , 
<  P2 


B. 


(7  -  1)U2' 


(3.47) 

(3.48) 

(3.49) 

(3.50) 

(3.51) 

(3.52) 

(3.53) 

(3.54) 


The  as  in  Duck  et  al.  are  to  be  determined  from  the  freestream  conditions.  In 

particular  in  order  to  write  the  exponential  terms,  arising  from  the  upstream  solution,  in 
terms  of  Bessel  functions  the  following  is  particularly  useful: 


exp[z(a;i  cot  /?  +  q;2)(cos  0  +  A  sin  0)1/2]  =  e  (3.55) 


where 


zL/  1  I  X  Ji'\J^e^y2h 


e  2 


z/=0 


1  +  Kfi 


^  .  _1  f  (<^1  +  “2)(cos0  +  Asin0)  + 


_  2 
2  ^*5 


ke^ 


(3.56) 


(3.57) 


Other  details  are  routine,  and  omitted  for  reasons  of  brevity. 

The  combined  system  (3.44)“(3.47)  then  represents  a  closed  (recursive-type)  system  de- 
termining  the  ai,+i, 
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4  Some  Numerical  Results 


The  first  set  of  data  we  present  is  for  the  particular  case  Mj  =  5,  0  =  25°  (weak  shock 
solution),  for  the  particular  case  of  a  slow  acoustic  mode  with  ax  =  1,  02  =  0,  and  a^  =  1. 
Fig.  la  shows  results  for  the  perturbation  pressure  on  the  wall  =  0).  Here,  and  in  all 
cases  the  solidus  denotes  the  real  part  of  a  function,  the  dashed  line  the  imaginary  part. 
This  figure  presents  a  picture  reminiscent  of  results  found  in  a  number  of  corresponding  two- 
dimensional  cases  by  Duck  et  al.  ^ ,  namely  that  of  an  oscillatory-type  nature  downstream. 
Fig.  lb  shows  the  corresponding  perturbation  shock  locationj  this  too  takes  on  a  similarly 
oscillatory  nature.  Figs.  2a,  2b  show  the  corresponding  results  for  the  fast  mode  case  (other 
data  remaining  the  same  as  for  figs.  1).  These  figures  are  to  be  compared  with  fig.  2  and 
fig.  3  of  Duck  et  al.  where  the  same  conditions  apply  except  that  a^  =  0  in  those  graphs. 
The  difference  between  the  two-dimensional  results  and  the  three-dimensional  results  are  in 
whether  or  not  the  disturbance  produces  an  acoustic  wave  behind  the  shock  that  intersects 
with  the  surface.  In  response  to  this  particular  slow  mode,  the  pressure  response  at  the 
surface  for  the  two-dimensional  disturbance  decays  algebraically  while  the  shock  position 
shows  a  characteristic  single  wavelength  response  in  the  far-field  limit.  On  the  other  hand, 
the  surface  pressure  response  to  the  three-dimensional  slow-mode  disturbance  does  not  decay 
and  shows  a  multi-wavelength  pattern  as  is  to  be  expected  if  a  decomposition  of  the  far  field 
behaviour  applies  to  the  three-dimensional  case  as  it  does  for  the  two-dimensional  case. 
The  shock  position  shows  a  two-wavelength  structure  as  would  be  expected  by  the  primary 
acoustic  disturbance  reflecting  off  the  wedge  surface  and  interacting  with  the  shock  from 
behind.  There  is  very  little  difference  between  the  responses  to  a  two-dimensional  fast¬ 
mode  disturbance  and  a  three-dimensional  fast-mode  disturbance.  A  search  in  parameter 
space  would  lead  to  the  same  conclusion  as  in  the  two-dimensional  case;  either  the  surface 
pressure  has  a  non-decaying  component  and  an  algebraically  decaying  component  or  it  has 
only  an  algebraically  decaying  component.  This  is  significant  in  that  the  solution  in  the 
absence  of  a  wedge  indicates  that  the  choices  are  between  only  a  non-decaying  response  and 
an  exponentially  decaying  response  (see  Hussaini  et  al.  1®,  Jackson  et  al.  Lasseigne  k 
Hussaini  This  issue  was  pursued  in  detail  in  Duck  et  al.  ^  where  it  is  determined  that 
the  condition  that  the  shock  remain  attached  to  the  apex  of  the  wedge  is  the  source  for 
the  algebraically  decaying  portion  of  the  response  -  the  apex  is  in  effect  a  point  source  for 
additional  acoustic  waves.  The  interaction  of  these  algebraically  decaying  acoustic  waves 
with  the  shock  in  turn  produces  an  algebraically  decaying  portion  in  both  the  entropy  and 
vorticity  modes.  Another  reason  that  it  is  important  to  determine  that  the  pressure  has 
only  a  decaying  or  sustained  oscillatory  behaviour  downstream  is  that,  as  mentioned  in  the 
analysis  of  the  response  to  the  two-dimensional  disturbances,  the  potential  functions  have 
components  of  exponential  growth  in  them.  The  calculation  of  the  physical  quantities  of 
pressure,  velocity  and  vorticity  are  therefore  dependent  on  some  fortuitous  cancellations,  or 
they  themselves  would  be  exponentially  growing.  In  the  following  section  we  investigate  the 
downstream  behaviour  of  the  three-dimensional  perturbation  solution  in  some  detail. 


5  The  Far  Downstream  Behaviour/Stability  Problem 


The  results  of  the  previous  section  show  that  a  general  growth  in  physical  quantities,  as  X2 
increases  does  not  occur  even  though  the  potential  functions  have  an  exponentially  growing 
component  in  them.  This  is  similar  to  the  two-dimensional  results  as  considered  by  Duck 
et  al.  in  which  it  was  shown  by  examining  the  recurrence  relations  that  only  decaying  or 
sustained  physical  solution  oscillations  were  possible.  The  same  type  of  analysis  is  possible 
for  the  three-dimensional  case,  but  indeed  has  its  own  peculiarities  as  will  be  seen. 

In  order  to  analyze  the  X2  oo  behaviour,  we  consider  the  homogeneous  solution  of  the 
system  (3.44)-(3.47),  and  consider  the  limit  dis  v  00.  In  particular,  we  write 

d^^DoK’',  and  a,  ^  AoK''e-''^\  (5.1) 


where  K  is  an  eigenvalue  whose  precise  value  (including  location  in  complex  space)  serves 
to  determine  the  ultimate  downstream  behaviour  of  the  disturbance.  If  |Ar|  >  1,  then  using 
the  generating  function  for  Bessel  functions,  we  have  (for  example) 


EZi  duMz)  ^  J:Z-ooDoK‘'Mz) 


(5-2) 

«  as  j^l  -4  00. 

Thus  exponential  growth  occurs  as  \z\  00  if  (aT  -  |  >  0.  If  we  substitute  (5.1) 

into  (3.44)  -  (3.47),  allow  u  00,  and  discard  the  inhomogeneous  terms,  then  we  obtain  the 
following  five  linear  homogeneous  equations  for  the  coefficients  Aq,  and  Do'- 


i.M  -  -  k 


+ k 

U2 


1;..  1. 

[  kU2p'^  ^  2^^^  k\ 


-  0{ui  -  U2)Do  =  0, 

af =  0, 


cos{^ 

+k 


2iu)U2 

+ 


4  [h/3^a,  K 


,-00  a: 


Ao 


iXcbUs  1;.  1 


+  sm(0  -  0)  +  e->‘K)Ao  +  +  ^cf>| 

-{u,  -  u,)  -  m,  {-g^  +  f  (/<  -  ^)}  Co  =  0. 

/  U  w/3  m  [  (2iu;U2  ,  e'’<>  „  \  ^ 


(5.3) 

(5.4) 


(5.5) 
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,  + Uix  -  4) 


COS 


(ff-6) 


M(^  +  e-»-A-M„  +  i/3CW  +  |cf> 


+ 


A2 

B2 


+ 


Vi  sin(/3  -  do) 


i^U2\  i 


Do 


A3 

B3 


^  2  \fcg/52  as  K  as  )  ] 


(5.6) 


This  equation  is  nonlinear  in  K  but  can  be  converted  to  a  tenth-order  linear,  generalized 
eigenvalue  problem,  which  may  be  solved  using  the  QZ  algorithm  by  using  the  additional 
variables 

io  =  KAo,  =  KCi^\ 

(5.7) 

Do  =  KDo. 

It  was  shown  by  Duck  et  al.  that  in  the  two-dimensional  case,  for  situations  in  which  the 
downstream  flow  was  supersonic,  there  were  four  imaginary  eigenvalues  (all  with  \K\  <  1), 
and  two  complex  eigenvalues  given  by 


(5.8) 


However,  although  these  complex  eigenvalues  are  such  that  \K\  >  1,  in  both  cases  Do  =  0, 
and  due  to  some  “fortuitous”  cancellations,  all  physical  quantities  remain  bounded  as  X2 
00  even  though  it  is  determined  that  the  potential  functions  themselves  grow  exponentially. 

In  the  three-dimensional  case,  we  may  expect,  ten  eigenvalues  instead  of  the  six  in  the 
two-dimensional  case.  However  it  turns  out  generally  that  there  are  just  eight  distinct 
eigenvalues  (two  of  the  eigenvalues  being  double  eigenvalues).  Generally  we  find:  (i)  four 
imaginary  eigenvalues,  (ii)  one  complex  conjugate  pair  of  double  eigenvalues  of  the  form 
±Kr  +  iKi,  and  (iii)  one  complex  conjugate  pair  of  eigenvalues  also  of  the  form  ±Kr  -f  iKi- 

For  9  =  25°,  Mi  =  5  ai  =  1  and  0:2  =  0,  the  four  imaginary  eigenvalues  have  magnitude 
less  than  unity  even  as  as  increases  from  zero.  The  eigensolutions  associated  with  these 
eigenvalues  do  not  contribute  to  any  growth  in  any  of  the  quantities  calculated.  Most  in¬ 
teresting,  however,  is  the  variation  of  these  four  eigenvalues  as  the  wedge  angle  is  increased. 
All  four  eigenvalues  approach  K  =  i  at  the  maximum  wedge  angle. 

The  (imaginary  part  of  the)  eigenvalues  described  by  (i)  are  shown  in  fig.  3a,  whilst  the 
real  and  imaginary  parts  of  the  eigenvalues  described  by  (ii)  and  (iii)  are  shown  in  figs  3b, 
3c  respectively.  Here  we  have  just  shown  the  eigenvalues  with  the  positive  value  of  Kr  (the 
other  eigenvalue,  corresponding  to  —K^  -H  iKi  may  obviously  be  simply  deduced).  We  next 
consider  the  effect  of  varying  9,  for  the  parameters  Mi  =  5,  ai  =  1,  aa  =  0,  and  03  =  1.  The 
variation  of  Ki  of  family  (i)  is  shown  in  fig.  4a,  and  the  variation  of  the  real  and  imaginary 
parts  of  eigenvalues  described  by  (ii)  and  (iii)  is  shown  in  figs.  4b  and  4c  respectively.  Again, 
just  the  positive  values  of  Kr  are  shown.  It  is  seen  that  the  double  eigenvalues  described  by 
(ii)  are  the  only  ones  with  |A|  >  1  and  therefore  lead  to  potential  exponential  growth.  It 
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turns  out  that  an  analytic  description  of  the  eigenvalues  described  by  (ii)  is  possible.  The 
result  is 


atk^ 


,1/2 


(5.9) 


It  is  also  possible  to  obtain  analytic  results  for  one  set  of  eigen-coefhcients  corresponding  to 
these  eigenvalues,  namely 


(5.10) 


C'«  = 


it/2  fl 


UJ 


ia^Ao  + 


iXka, 


U. 


2  V  K 


Dq  —  0. 


(5.11) 

(5.12) 

(5.13) 


The  above  solution  readily  reduces  to  the  two-dimensional  solution  in  the  limit  cva  — )•  0.  Some 
understanding  of  the  origin  of  this  other  family  of  coefficients  for  this  eigenvalue,  and  also 
of  the  additional  complex  eigenvalue  (iii)  may  be  made  by  considering  the  two-dimensional 
limit,  03  -J-  0  and  seeking  a  solution  to  the  system  with  Aq  =  Do  =  =  0.  After  some 

algebra,  we  find  the  following  (four)  eigenvalues: 


hCU2 


Xikus 


-.2 


u 


keCU2 


1}^ 


(5.14) 


together  with 


U2kii 

2uj 


¥ 


2Xikas 

'uJI 


(5.15) 


Two  of  the  above  values  for  K  correspond  to  (5.9),  and  this  partly  explains  the  origin  of  the 
double  eigenvalues.  As  03  increases  from  zero,  it  is  expected  that  Aq,  Dq,  and  will  no 
longer  be  zero  for  this  branch.  The  other  pair  of  eigenvalues  correspond  to  the  eigenvalues 
described  by  (iii). 

Again  it  is  important  to  note  is  that  with  the  exception  of  the  eigenvalues  (ii),  all  eigenval¬ 
ues  of  K  correspond  to  |/T|  <  1,  and  thus  are  of  limited  physical  significance.  Interestingly  at 
the  maximum  value  of  9,  (corresponding  to  the  maximum  wedge  angle  for  attached  flow  at  a 
particular  freestream  Mach  number),  all  eigenvalues  take  on  the  same  value,  K  =  i.  Turning 
now  to  the  effects  of  the  eigenvalues  described  by  (ii),  again,  as  in  the  two-dimensional  case, 
it  turns  out  that  although  seemingly  these  eigensolutions  correspond  to  unstable  (grow¬ 
ing)  downstream  solutions,  in  fact  the  evaluation  of  physical  quantities  leads  to  solutions 
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downstream  which  remain  bounded.  This  is  best  illustrated  by  defining  the  following  set  of 
functions: 


^  =  AqC 


f:  K-Mkim), 

l/zz~00 

(5.16) 

oo 

(5.17) 

where  /T  is  defined  by  (5.9)  and  we  write  E  =  It  would  be  expected  that 

E  and  ^  determine  the  far-downstream  behaviour  of  the  flow  which  grows  downstream. 
However  inspection  of  the  solution  above,  in  the  cases  considered,  leads  to  the  conclusion 
that 

V(^=-VAE.  (5.18) 

This  then  leads  to  all  velocity  components  (and  hence  all  physical  quantities)  remaining 
bounded  downstream.  Numerically,  it  is  determined  that  all  physical  quantities  remained 
bounded  no  matter  what  the  parameters  selected  for  the  calculation,  thus  the  above  must 
hold  for  all  eigensolutions  with  \K\  >  1. 

In  summary,  therefore,  we  see  that  provided  the  flow  downstream  is  supersonic  then 
disturbances  will  not  grow  downstream.  This  boundary  of  downstream  behaviours  is  inde¬ 
pendent  of  the  crossflow  wavenumber  a^,  and  occurs  along  the  ’line’  where  ^  =  0  (where 
k  and  simultaneously  become  singular,  and  also  change  from  being  real  to  imaginary 
quantities).  The  result  is  that  all  of  the  Bessel  functions  in  the  series  (3.21)-(3.23)  will  grow 
downstream  if  the  downstream  flow  is  subsonic  (although  we  note  the  comments  of  Salas  k 
Morgan^ ^  stating  that  strong  shock  solutions  could  perhaps  be  stable  under  a  different  set 
of  boundary  conditions). 

We  now  address  the  nature  of  the  downstream  response  to  sustained  excitation  (i.e. 
acoustic  waves,  entropy  waves,  or  vorticity  waves).  Arguments  similar  to  those  employed  by 
Duck  et  al.  ^  are  again  appropriate.  The  frequency  of  the  disturbance  downstream  of  the 
shock  is  determined  by  the  frequency  of  the  disturbance  upstream  of  the  shock  and  both 
wavenumbers  of  the  disturbance  tangential  to  the  shock  are  also  fixed  by  the  deflection  of  the 
shock  produced  by  the  upstream  disturbance.  There  is  a  finite  range  of  tangential  wavenum¬ 
bers  that  exclude  the  existence  of  plane  acoustic  waves  which  propagate  downstream.  If 
the  tangential  wavenumbers  lie  outside  of  this  range,  then  a  plane  wave  propagates  at  an 
angle  Op  to  the  shock.  In  the  case  of  the  tangential  wavenumbers  being  within  the  excluded 
range,  an  acoustic  field  is  generated  that  decays  algebraically  rather  than  exponentially  as 
would  be  the  case  in  the  absence  of  the  wedge.  If  the  tangential  wavenumber  is  outside  of 
the  excluded  range,  then  the  solution  in  the  presence  of  the  wedge  exhibits  two  types  of  be¬ 
haviour,  dependent  upon  the  angle  of  propagation  of  an  acoustic  disturbance  in  the  absence 
of  the  wedge.  If  the  angle  of  propagation  Op  is  greater  than  the  angle  between  the  shock  and 
the  wedge  —  0,  then  the  pressure  disturbances  generated  at  the  shock  intersect  with  and 
reflect  from  the  wedge  surface,  leading  to  a  non-decaying  pressure  field  everywhere  between 
the  shock  and  the  wedge  surface,  with  a  superimposed  algebraically  decaying  pressure  field 
owing  to  the  requirement  of  an  attached  shock  and  zero  velocity  normal  to  the  wedge  sur¬ 
face.  If  the  angle  of  propagation  Op  is  less  than  the  angle  between  the  shock  and  the  wedge 
surface,  then  the  pressure  field  is  divided  into  two  regions  by  a  ray  emanating  from  the  apex 
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and  parallel  to  the  direction  of  propagation  of  the  pressure  disturbance  in  the  absence  of 
the  wedge.  Between  the  shock  and  this  ray  there  is  a  non-decaying  pressure  field,  with  the 
aforementioned  superimposed  algebraically  decaying  component.  Between  the  ray  and  the 
wedge  surface  the  pressure  field  is  algebraically  decaying  and  there  is  no  component  with 
sustained  oscillations.  This  was  referred  to  as  the  shadow  region  by  Duck  et  al.  ^ 

The  various  downstream  limits  are  determined  by  the  value  of  B.  If  6  is  real,  then  the 
wavenumber  of  the  imposed  disturbances  lies  within  the  excluded  range,  and  the  pressure 
along  any  ray  62  =  constant  decays  algebraically  For  62  =  Bq  the  solutions  does  not  decay 
leading  to  the  appropriate  oscillatory  behaviour  necessary  to  match  conditions  at  the  shock. 
If  B  is  complex,  then  the  condition 


(ai  cot  /3  -f  q;2)(cos  0  -|-  A  sin  0)  -f 

ki 


>  1, 


(5.19) 


is  satisfied.  We  also  note  that  for  equation  (3.55)  to  be  satisfied  then  sin  B  must  remain  real, 
which  requires  that 

B  =  ±it/2  -  iBi,  (5.20) 

where 

(  (cKi  cot  ^  +  02)  (cos  ^  -h  A  sin  ^)  -H  ) 

Bi  =  cosh-'  I  ± - \ .  (5.21) 

Bi  is  related  to  the  angle  between  the  direction  of  the  shock  and  the  direction  of  the  acoustic 
disturbance  Bp  through  tanh^,-  =  /3tan0p.  If  k  is  greater  than  ^o,  then  the  generated  or 
transmitted  acoustic  waves  intersect  the  wedge  surface.  For  0  <  <  Bq,  there  are  two 

regions:  a  region  of  sustained  oscillatory  behaviour  of  the  pressure  field  near  the  shock  when 
B2  >  Bi,  and  a  region  of  algebraic  decay  of  the  pressure  field  near  the  wedge  surface  when 

02<k 

These  features  are  common  to  those  found  in  the  two-dimensional  case  by  Duck  et  al. 
'.  Additionally,  we  can  expect  that  the  downstream  (forced)  behaviour  will  take  the  same 
form  as  in  the  two  dimensional  case.  Taking  the  wall  pressure,  for  example,  we  may  write 
the  decomposition 


(5.22) 

i.e.,  a  solution  with  an  oscillatory  component,  two  decaying  acoustic  waves,  and  a  (faster) 
decaying  component,  Q{x2). 


6  Conclusions 

The  interaction  between  a  shock  wave  and  three-dimensional  freestream  disturbances  has 
been  considered,  in  particular  the  nature  of  the  flow  far  downstream.  The  somewhat  fortu¬ 
itous  cancelations  that  were  found  to  occur  in  the  analogous  two-dimensional  work  (Duck 
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et  al.  again  intriguingly  occur,  yielding  a  non-growing  physical  solution  downstream, 
thus  confirming  the  stability  of  shocks  with  downstream  supersonic  flow,  in  line  with  the 
widely-held  belief  regarding  the  stability  of  such  shock  waves.  Our  results  also  point  to  the 
‘instability’  of  shocks  with  downstream  subsonic  flow,  since  in  that  case  the  argument  of 
the  Bessel  functions  in  (3.21)-(3.24)  becomes  imaginary,  and  as  a  result  exponential  growth 
will  be  expected  downstream.  Finally  we  note  that  the  location  of  the  boundary  between 
growing  and  non-growing  downstream  waves  is  unaffected  by  the  crossflow  wavenumber  03. 
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Figure  2a.  Wall  pressure  perturbations,  Mi  =  5,  0  =  25°,  (weak-shock  solution),  fast 

acoustic  mode,  =  l,a2  =  0,a3  =  1. 
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